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An application of mnltivariate total positivity to 

peacocks 

ANTOINE MARIE BOGSO * 


Abstract: We use multivariate total positivity theory to exhibit new families of peacocks. As 

the authors of |18], our guiding example is the result of Carr-Ewald-Xiao |9]. We shall introduce the 
notion of strong conditional monotonicity. This concept is strictly more restrictive than the conditional 
monotonicity as defined in [18] (see also [2], [5] and [36]). There are many random vectors which are 
strongly conditionally monotone (SCM). Indeed, we shall prove that multivariate totally positive of order 
2 (MTP 2 ) random vectors are SCM. As a consequence, stochastic processes with MTP 2 finite-dimensional 
marginals are SCM. This family includes processes with independent and log-concave increments, and 
one-dimensional diffusions which have absolutely continuous transition kernels. 

Key words: convex order, peacocks, total positivity of order 2 (TP 2 ), multivariate total positivity of 
order 2 (MTP 2 ), Markov property, strong conditional monotonicity. 


1 Introduction 

We call peacock a real valued process {Zt,t > 0) which is iiitegrable, i.e. ]E[|^t|] < oo for 
every t > 0, and which is increasing in the convex order, i.e., for every convex fonction 
'(/i : R —>■ K, 


the map: t G M+ '—> E[ijj{Zt)] G] — oo, +oo] is non-decreasing. (1-1) 

Observe that if {Zt, t > 0) is a peacock, then K[Zt] does not depend on t. Indeed, it suffices 
to apply (1.1) first with ^/’(a;) = x, then with The pun peacock comes from 

the french: “Processus Croissant pour I’Ordre Convexe" which acronym: “P.C.O.C" may 
be pronounced “peacock". To prove that an integrable process satisfies (1.1), it suffices to 
consider convex functions which belong to the set: 

C := {'0 : R —t R convex C^-function such that ip" has a compact support}, 

where 0" denotes the second order derivative of 0. Note that if 0 G C, then its derivative 
0' is bounded and there exist fci, A ;2 > 0 such that: 

VxGR, |0(a:)| < fci-I-A:2|a:|. 

There are two remarkable results which motivate the investigation on peacocks. In 2008, 
Carr-Ewald-Xiao [9] proved that Asian options are increasing in the convex order with 
respect to the maturity. In other words, if {Bs,s > 0) denotes a standard Brownian 
motion issued from 0, then 

:= - y is a peacock. (CEX08) 
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The second result due to Kellerer [28] states that a real valued process {Zt,t > 0) is a 
peacock if and only if there exists a martingale {Mt ,t > 0) with the same one-dimensional 

marginals as {Zt,t > 0), i.e., for every t > 0, Mt Zt- This martingale is not unique 
in general and it may be chosen Markovian. Recently, Hisch-Roynette [19] offered a new 
proof of Kellerer’s theorem. However, the Kellerer’s proof is not constructive, and then 
it helps neither establishing whether or not a process is a peacock, nor constructing an 
associated martingale to a given peacock {Zt, t > 0), i.e. a martingale having the same one¬ 
dimensional marginals as {Zt,t >0). In [7], Baker-Yor provide an associated martingale 
to (Nt,t > 0) using the Brownian sheet. Inspired by Carr-Ewald-Xiao and Baker-Yor 
results, the authors of [18] exhibited several examples of peacocks and they provided 
several methods to associate explicitely martingales to certain of them. We refer the 
reader to [32], [17], [3], [8] and [6] for further interesting results about peacock processes. 
In this paper, we exibit new families of peacocks using multivariate total positivity theory. 
But for many of them, finding an associated martingale remains open. Let us mention 
that total positivity is a nice property that plays an important role in various domains 
of Mechanics and Mathematics. There is a large amount of literature concerning total 
positivity. We shall follow Karlin [24] and Karlin-Rinot [29] for basic definitions and 
results. 

In section 2, we give some basic results concerning total positivity and multivariate total 
positivity of order 2. Section 3 is reserved to strong conditional monotonicity results. 
Finally, in section 4, we use strong conditional monotonicity theorems to exhibit new 
classes of peacocks, inspired from the Carr-Ewald-Xiao example. 


2 Total positivity and multivariate total positivity 
of order 2 

We first define totally positive functions of order 2 and give several examples of Markov 
processes with totally positive transition kernels. Then, we deal with an extension of total 
positivity of order 2 to multivariate distributions. 


2.1 Totally positive functions of order 2 

We follow the terminology and notation of Karlin [23]. 

Definition 2.1. A function p : R x K. —>■ R.+ is said to be totally positive of order 2 (TP 2 ) 
if for every real numbers xi < X 2 , yi < 2 / 2 , 


P 


f Xi,X2 \ 

:= det ( 

V 2/1,1/2 ) 

\ 


p{xi,yi) p{xi,y 2 ) 


> 0 . 


(TP2) 


pix 2 ,yi) p{x 2 ,y 2 ) 

Similarly, a function p : Z x Z — >■ R+ is said to be TP2 if, for every integers ki < k2 and 
h < h, P satisfies (TP2). 

Note that one may define totally positive functions of order higher than 2 (see [23]). 
Remark 2.2. Let D denote a subset of R x R. which satisfies the following property: 

For every xi < X 2 and yi < y 2 , I 

\ (P) 

[{xi,y2) G D and {x2,yi) G D] [{xi,yi) G D and (x 2 ,p 2 ) G L>]. J 

Let p : D ^ R+ be TP2, i.e. for every xi < X2, yi < P2 such that (xi,yi), (xi,y2), (x2,yi) 
and {x2,y2) belong to D, 


Xl,X2 

2 / 1 , 2/2 


> 0 . 
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We define: 


Vx, j/ G R, p{x,y) 


p{x,y) if{x,y)eD, 
0 otherwise. 


Then, p is TP 2 if and only if p is TP 2 . 

Here are some examples of Z? C R x R satisfying (P). 

i) If I and J are two intervals of R, then I x J satishes (P). 

ii) For every reals ko < ki and every (a, /3) G R x R \ {(0, 0)}, 

D = {{x,y) & ko < ax — jSy < ki} satifies {P). 


We give properties of TP 2 functions assuming that they are defined on R x R. By the 
preceding remark, one may extend these results to functions defined on subsets of R x R 
which satisfy (P). 

The following characterization result of smooth TP 2 functions is proved in [24]. 


Proposition 2.3. (Karlin [24])- Let p : R x R — >■ R 

dp dp , d'^p , / ^ f m, 

—, —, and exist at each point {x,y) 0 / R x 

dx dy dxdy 

1) If p is TP 2 , then, for every reals xi < X 2 and y, 


he such that the partial derivatives 


det 


( Pixi,y) \ 

p{x2,y) ^{x2,y) 

dy 


> 0 , 


( 2 . 1 ) 


and, for every (x, j/) G R x 


det 


p(x,y) 






> 0 . 


( 2 . 2 ) 


2J Conversely, if p{x,y) > 0 for every {x,y) G R x R, then (2.2) implies (2.1), which in 
turn implies that p is TP 2 . 

A second criterion of smooth TP 2 functions follows from Proposition 2.3. 


Corollary 2.4. (Karlin [24]). Let p : R x 
each point {x, y) o/R x 
p is TP 2 if and only if 


->■ R+ be strictly positive and such that at 
d^ (logp) 

each point {x,y) o/R x R, the second order partial derivative —— iK^y) exists. Then 


dxdy 


d^jlogp) 

dxdy 


{x,y) > 0. 


Proof. This follows from (2.2) and from the straightforward relation: 
,92 (logp) 


P {x,y)- 


dxdy 


-{x, y) = det 


f p{x,y) \ 




□ 
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Example 2.5. (Brownian transition densities). 
Consider the family (pt : K. x K. —>■ R+, t > 0) given by: 


For every t > 0, 


V (x, j/) G R X R, pt{x, y) = exp 


(92(logpt) 1 


-{x - yf 

2t 


and, by Corollary 2.4, pt is TP 2 . 

More generally, if / : R —>■ R+ is a strictly positive C^-function, then (x, y) 1 —/(x — y) is 
TP 2 if and only if / is log-concave. Indeed, 


a" 

dxdy 


[log /(x - y)] = -(log /)"(x - y). 


Example 2.6. (Ornstein-Uhlenbeck transiton densities). 
Let {pt,t > 0) be the densities defined on R x Rby: 


Ptix,y) = 

For every t > 0, 


-——^—r exp —ce 
27r sinh(d:) \ 

d'^ {log Pt) 


; (j/- xe-=* - 1^(1 - 6-^=*))^ 
2 sinh(ct) 


(c, ly G 


{x,y) = 


> 0 , 


dxdy sinh(ct) 

and, from Corollary 2.4, pt is TP 2 . 

Counterexample 2.7. The function p defined by: 

Vx, 2 /GR, p{x,y) = ^ -— 

1 -h {x-y)^ 

is not TP 2 . 

New TP 2 functions are generated using the following classical composition formula: 

be two Borel functions, and let a be a positive measure 


Lemma 2.8. Let p, q : R x 
on R such that: 


Vx, z G 


\p{x,y)\\q{y,z)\ a{dy) < 00 . 


Let r denote the function defined on R x R by: 

Vx, 2 GIR, r{x,z)= / p{x,y)q{y,z)a{dy). 

Jr 

Then, for every reals xi < X 2 , zi < Z 2 , 


Xi,X2 

21,^2 


P 


yi< 


Xi,X2 


<y 2 V 


q 


yi,y2 

Zl,Z2 


cr{dyi)a{dy 2 ). 


(2.3) 


Remark 2.9. If p and q are two integrable functions with respect to the Lebesgue measure 
and li r := p * q denotes the convolution product of p and q, i.e. 


VxGR, r(x) = / p{x - y)q{y)dy, 

Jr 
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then, using (2.3), for every xi < X 2 and zi < 22 , 
r{xi - zi) r{xi - Z 2 ) 


det 


r{x2 - zi) r{x2 - Z 2 ) 

p{xi - yi) p{xi - 2/2) 


det 


det 


Vl<V2 


p{x2 - yi) p{x2 - 2/2) 


9(2/1 - zi) 9(2/1 - Z2) 
9(2/2 - zi) q{y2 - Z2) 


dyidy2. 


From Lemma 2.8 and Remark 2.9, we easily deduce the following result: 

X M —>■ ]R_|_ be two TP 2 functions such that the product 


Proposition 2.10. Let p^q : 

r : R X R —>■ R_|_ given by: 


r{x,z)= / p{x,y)q{y,z)dy 

Jr 


is finite. Then r is TP 2 . 

In particular, if p and q are two integrable log-concave functions, then the convolution 
product r = p * q is also log-concave. 

Remark 2.11. Proposition 2.10 allows to regularise TP 2 functions in such a way that 
total positivity is preserved. Indeed, if g is a TP 2 function such that, for every z G R, 
q{-,z) is integrable, then, for every e > 0, the function q^ defined by: 


Va;, z G R, qe{x, z) 



[ 

{x-yf 

/ exp 
/B 

[ 2e2 J 


qiy,z)dy 


is TP 2 . Moreover, 


VzGR, lim g£(-, z) = g(-, z) in L^(R). 


2.2 Markov processes with totally positive transition kernels 

In this paragraph, we present some examples of Markov processes with totally positive 
transition kernels. 

Definition 2.12. Let P := {Ps^t{x,dy),Q < s < t, x G I) be the transition function of a 
Markov process {{Xt,t > 0), {Px,x G I)) with values in a sub-interval I o/R. P is said to 
be totally positive of order 2 (TP 2 ) if, for every 0 < s < t, every xi < X 2 elements of I, 
and every Borel subsets Ei, E 2 of I such that Ei < E 2 (i.e. oi < 02 for every oi G Ei 
and 02 G E 2 ), we have: 


P. 


Pf 


Xl,X2 

El, E 2 


Xl,X2 

El, E 2 


:= det 


> 0 . 


Ps^t{xi,Ei) Ps,t{xi,E2) 

Ps,t{x 2 ,Ei) Ps,t{x2,E2) 

Suppose moreover that {Xt,t > 0) is time-homogeneous. Then P is TP 2 if and only if 

Pt{xi,Ei) Pt{xi,E 2 ) 

Pt{x2,Ei) Pt{x2,E2) 


(2.4) 


:= det 


> 0 . 


(2.5) 


Remark 2.13. Let P := {Ps^t{x,dy),0 < s < t,x G I) denote the transition function 
of a Markov process {{Xt,t > 0), {Vx,x G /)) taking values in a sub-interval / of R. We 
suppose that, for every 0 < s < t and x G I, Ps,t{x, dy) has a continuous density Ps,t{x, •) 


5 



ANTOINE MARIE BOGSO 


with respect to a cr-finite regular measure. Then P is TP 2 if and only if is TP 2 , i.e. 
for every xi < X 2 , yi < 2/2 elements of I, 


Ps,t 


f Xi,X2 \ 

:= det ( 

V VTV2 ) 

\ 


Ps,t{xi,y2) 
Ps,t{^2,yi) Ps,ti^2,y2) 


> 0 . 


Definition 2.14. Let P := {Pg^tik, l),0 < s < t, (fc, 1) G I x I) be the transition function 
of a continuous time Markov chain {Xt, t > 0) which takes values in a sub-interval I ofL. 
We say that P is TP 2 if, for every 0 < s < t and every integers ki < k 2 , h < I 2 in I, 




ki,k2 

ll,l2 


:= det 


Ps,t{ki,li) Psyikijh) 
Ps,t{k2,h) Ps,t{k2,l2) 
If (Xt,t > 0) is time-homogeneous, then (2.6) is equivalent to: 

PtikiJi) Ptiki,l2) 


> 0 . 


( 2 . 6 ) 


Pt 


ki,k2 

h, I 2 


:= det 


> 0 . 


Pt{k2,h) Pt{k2,l2) 


(2.7) 


There are many Markov processes with totally positive transition kernels. Let us give 
some of them. 


2.2.1 Processes with independent and log-concave increments 

Let {Xt,t > 0) be a real valued process with independent increments, i.e. 

VO < s < t, Xt — ATs is independent of Fg '■= o'(X„; 0 < u < s). (PH) 

We suppose that, for every 0 < s < t, the increment Xt — Xg is log-concave. In other 
words, Xt — Xg has a density (with respect to Lebesgue measure) which is log-concave, 

i.e. 

yx,y G K, 9 g]0, 1[, pg^tiOx -h (1 - 9)y) > (ps,t(a;))®(ps,t( 2 /))^"^. (2.8) 

If Xt — Xg takes values in Z, then (2.8) may be replaced by 

pl,t{k) > Ps,t{k - l)pg,t{k + I), (2.9) 

where, for every fc G Z, pg,t{k) := P(Xt — Xg = k). 

Many common r.v.’s are log-concave. Indeed, Gaussian, uniform, exponential, binomial, 
negative binomial, geometric and Poisson r.v.’s are log-concave. On the contrary. Gamma 
r.v.’s with parameter a G (0,1) are not log-concave. 

Theorem 2.15. (An [1], Daduna-Szekli[10]). A Lebesgue-measurable function / : K. —>■ 
1R+ is log-concave if and only if 

( fixi - yi) f{xi - 2 / 2 ) \ 

V xi < a ;2 G M, 2/1 < 2/2 G R, det >0. (2.10) 

\ f{x2 - 2/1) f{x2 - 2/2) / 

Then, {Xt,t > 0) is a Markov process with transition function Pg^t{x,dy) given by: 

V0<s<t, xGR, Pg,t{x, dy) = pg,t{y - x)dy, (2.11) 

and, since Xt — Xg is log-concave, we deduce from (2.11) and from Theorem 2.15 that 
Pg,t{x,dy) is TP 2 . 
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2.2.2 Absolute value of a process with independent, symmetric and PFqo 
increments 

Let {Xt,t > 0) be a real valued process with independent increments (PII) such that, for 
every 0 < s < t, Xt — Xg is symmetric and has a density denoted ps^t- Note that Ps^t 
is symmetric, i.e., for every x G K+, Ps,t{x) = Ps,t{—x). Then, (|Ari|,t > 0) is a Markov 
process whose transition function is given by: 

dy) := y)dy (0 < s < t, a; G M), 

where 

p*sM^ y) = -x)+ PsA-^ - y)- 

We suppose in addition that, for every 0 < s < t, Ps,t is a Polya frequency (PFoo) 
function, i.e. for every integers r > 1, 1 < m < r, and every a;i < 0:2 < • • • < Xm, 

yi < y 2 < ■ ■ ■ < y-m, 

det[psAxi - yj); 1 < *,J < to] > 0. (PFoo) 

Then the transition function PsAx, dy) is TP 2 , i.e., for every 0 < s < t, p* ^ is TP 2 . This 
is a direct consequence of the following result due to Karlin [23]. 

Theorem 2.16. (Karlin [23]). Let / : M —>■ 1R+ be a symmetric PFoo density function. 
Then the function f* defined on R_|_ x R_|_ by 

f*{x, y) = fix -y)+ fi-x - y) 


is a TP 2 function. 

We mention that Theorem 2.16 remains valid if we consider discrete and symmetric 
PFoo densities. To prove Theorem 2.16, the author applies Schoenberg’s characterisation 
of symmetric PFoo densities in terms of their Laplace transforms. 

Theorem 2.17. (Schoenberg [35]). A symmetric density function / : M —>■ M+ is PFoo 

/ OO 

fiy)dy exists in a strip (of the 

-OO 

complex plane) including the imaginary axis in its interior and has the form 


$(s) = 


n(l-afs2) 

2=1 


( 2 . 12 ) 


where a > 0, G M for every i, and 0<a+'^ a] <co. 

i=l 

A discrete analog of Theorem 2.17 has been proved by Edrei [11]. 

Theorem 2.18. (Edrei [11]). A symmetric density function f : 1, ^ R+ is PFoo if and 

OO 

only if its Laurent series S{z) = /(fc)z* converges in some ring (of the complex 

k——oo 

plane) including the unit circle in its interior, and the analytie continuation S of S is of 
the form 

n ((l + af)+a.(z+i)) 

^(z) = Cexp (a (z + z ^)) LJ—-1, (2.13) 

n(l-7i2)(l-7*2-1) 

2=1 

OO 

where C > 0 , a > 0 , > 0 , 0 < 7 ^ < 1 for every i, and Y) Ai + 7 *) < 00 . 

i=l 
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Here are some examples of symmetric PFqo density functions. 


Example 2.19. By Theorem 2.17, 

i) a symmetric gaussian density is PFqo, 

ii) for every A > 0, the density function f{x) = ^ exp(—A|a;|) is PFoo- 


Example 2.20. Theorem 2.18 applies in the following cases: 
i) For every a G M+, the density function / : Z —>■ 1R+ defined by 

1 + q:^ 


/( 0 ) = 


/(I) = /(-!) = 


a 


(l + a)2’ (l + a)2’ 

and f{k) = 0 if fc ^ {~1) 0; 1} is PFoo- Indeed, if we define 


OO 

S{z):= ^ 

k— — oo 


then 


S{z) = 


(l + a)2 


(1 + a^) + a I z + 


which is of the form (2.13). 

(1 — 6*^)2 

ii) Let a > 0 and c =-The density function / : Z —>■ ]R_|_ given by 

1 — e“2a 

/(O) = c and f{k) = fc G Z \ {0} 

OO 

is PFoo, since its Laurent series ^ f{k)z^ admits the representation 

k——oo 


S{z) = 


(l-e-“)2 


(1 - e-“0) 1 - 


which is of the form (2.13). 

2.2.3 One-dimensional diffusions 

An important class of real valued Markov processes with TP 2 transition kernels consists 
of one-dimensional diffusions. 

Theorem 2.21. (Karlin-Taylor [27], Chapter 15, Problem 21). 

Let {{Xt,t > 0), (Pa;,x G I)) be a one-dimensional diffusion on a sub-interval I o/R, and 
let {Pt(x,dy),t > 0,x G I) be its transition function. We suppose that (t,x) 1 —> Pt{x,dy) 
is continuous in x for every t. Then, for every t > Q, Pt{x,dy) is TP 2 , i.e. for every 
xi,X 2 G I and every Borel subsets Ei < E 2 of I, 

/ Pt{xi,Ei) Pt{xi,E2) \ 

det > 0. (2.14) 

\ Pt{x2,Ei) Pt{x2,E2) ) 

In particular, if Pt(x,dy) = pt(x,y)dy, and if pt is continuous in y for every x, then pt is 
a TP 2 function. 
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This result is due at origin to Karlin-McGregor [25] who showed that (2.14) has a prob¬ 
abilistic interpretation: Suppose that two particles qi and q 2 , started at time zero in states 
xi and X 2 respectively, execute the process {Xt,t > 0) simultaneously and independently. 
Then the determinant 

/ Pt{xi,Ei) Pt{xi,E2) \ 

det (2.15) 

\ Pt{x2,Ei) Pt{x2,E2) ) 

is equal to the probability that at time t, qi is located in Ei and q 2 is located in E 2 without 
these particles having occupied simultaneously a common state at some earlier time t < t. 
Karkin-McGregor [25] reach the same result for several time-homogeneous and strong 
Markov processes whose state space is a subset of the real line. In particular, the transition 
probability matrix of a birth and death process is TP 2 . 

Moreover, as a consequence of Theorem 2.21, bridges of one-dimensional diffusions have 
TP 2 transition functions. We refer to [15] for a rigorous definition of the bridge of a 
one-dimensional diffusion. 

2.3 Multivariate total positivity of order 2 

Definition 2.22. 

1) A function p : M” —>■ M-i- (n G N, n > 2^ is said to be multivariate totally positive of 
order 2 (MTP 2 ) if for every x = {xi, • • • , Xn) and y = (j/i, • • • , ?/„) in R", 

p(x A y)p(x V y) > p(x)p(y), (MTP 2 ) 

where 

xAy = (min(a;i,?/i),--- ,min(a;„,y„)) 
and (2.16) 

xVy = (max(a;i,yi), • • • , max(a:„, ?/„)). 

2) A random vector {Xi,--- ,Xn) with real components is said to be multivariate totally 
positive of order 2 (MTP 2 ) if it is absolutely continuous with respect to a a-finite product 
measure (which we shall always denote dxi ■ ■ ■ dxn) and if its density p : R" —>■ R_|_ is 
MTP 2 . 

Remark 2.23. By definition, the MTP 2 property is invariant under permutations, i.e. 
if a random vector (ATi, • • • , X„) is MTP 2 , then, for every permutation tt of {1, • • • , n}, 
(-’^7r(l), • • • W7i-(n)) is MTP 2 . 

The following result is proved in Karlin-Rinot [29]. 

Theorem 2.24. (Karlin-Rinot [29]). 

1) If {Xi,--- ,Xn) is a MTP 2 random vector, then, for every k G N, 2 < k < n, 
[X,,--- ,Xk) IS MTP 2 . 

2) Let p : R” —>■ R+ be a MTP 2 density. Then, for every k G^, 2 < k < n, and for every 
continuous and bounded functions /^ : R —>■ R_|_, i = 1, • • • ,n, the function p^^l : R^ —>■ R_|_ 
given by: 

P^'^'>{xi,- ■■ ,Xk) 

k p n 

= T\Mxi) p{xi,-■ ■ ,Xk,Uk+l,-■ ■ ,Un) TT fj{Uj)dUk+l---dUn 

^=l j=k+l 

is MTP 2 . 

As a consequence of Theorem 2.24, we have: 
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Corollary 2.25. (Karlin-Rinot [29]). 

Let m, n G N*. Let p -.'S) x R™ —>■ R_|_ and q : R™ x R" —>■ R_|_ be two MTP 2 densities. 
Then the function r : R* x R" defined by: 

V(a;,z) G R' X R”, r(a;,z) = [ p{x,y)q{y, z)dy 


is MTP 2 . 

Remark 2.26. Corollary 2.25 allows to regularize MTP 2 densities while preserving MTP 2 
property. 

Note that MTP 2 random vectors satisfy a co-monotony principle. Indeed, Sarkar [34] 
proved the following result. 

Theorem 2.27. (Sarkar [34]). Let X := {Xi,--- ,X„) (n G be MTP 2 , and let 
(j),ip : R” —>■ R &e two measurable and simultaneously (componentwise) non-decreasing 
(resp. non-increasing) onR". Then 

E[(/.(X)(^(X)] > E[</>(X)]E[(^(X)]. (2.17) 

Recently, Pages [31] introduced a functional co-monotony principle for stochastic pro¬ 
cesses X := {Xt,t > 0). He offered an extension of Inequality (2.17) by replacing <j) and 
ip with functionals of the hole path of X. The author also provided many examples of 
co-monotone stochastic processes. 

Here are some basic examples of MTP 2 distributions (see Karlin-Rinot [29]). 

Example 2.28. (Karlin-Rinot [29], Section 3). 

1) If Xi, • • • , Xn is a sample of i.i.d. random variables, each Xi having a density, then the 
joint density of the order statistics • • • , is MTP 2 . 

2) A Gaussian random vector (Xi, • • • , X„) with an invertible covariance matrix S is MTP 2 
if and only if the inverse matrix of E has negative off-diagonal elements. 

3) Let (Xi, • • • , Xn) be a gaussian random vector with zero mean, and with an invertible 
covariance matrix E. Let E“^ denote the inverse matrix of E. Then (jXij, • • • , |X„|) is 
MTP 2 if and only if there exists a diagonal matrix D with elements ±1 such that DT,~^D 
has negative off-diagonal elements. 

4) Let {Xt,t > 0) be a Markov process with absolutely continuous and TP 2 transition 
kernel. Then, for every distinct elements ti, • • • ,tn in R+, (Xj^, • • • , Xt„) is MTP 2 . 

Further examples of MTP2 distributions may be found in Karlin-Rinot [29] and Gupta- 
Richards [16]. 


3 Strong conditional monotonicity 

We introduce the notion of strong conditional monotonicity which strictly implies condi¬ 
tional monotonicity as defined in [5] and [18]. 

Definition 3.1. For every n G N*, let denotes the set of continuous and bounded 
functions (j) : R” —>■ R which are componentwise non-decreasing. 

Definition 3.2. (Strong conditional monotonicity). 
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An application of multivariate total positivity to peacocks 


1) A random vector (Xi,--- ,Xn) with real components is said to be strongly condition¬ 
ally monotone (SCM) if, for every i G {I,-"’ every continuous and strictly positive 
functions /^ : M —>■ M+, k = 1, • • • ,n such that: 


and every (j) G In, we have: 


z G 


Ki{n,z) 


E 


n 


< 00 , 


E 

1_1 

n 

•• n fkiXf,) 

k=l 

1 - 1 

II 

E 

1_1 

X, =z 



(3.1) 


(SCM) 


is a non-decreasing function. 


2) A real valued process {X\,X > 0) is said to be strongly conditionally monotone (SCM) 
if its finite-dimensional marginals are SCM. 

Remark 3.3. 

1) If there are subsets /i,--- ,In of R such that, for every k, Xk takes values in Ik, we 
may suppose in Definition 3.2 that (j) is defined on /i x • • • x Ik, fk is defined on Ik, and 
2 I— > Ki{n, z) is dehned on A. 

2) Let (Xa, a > 0) be a real valued process, and let 0 : IR+ x R —>■ R be such that, for every 
A > 0, a; I—^ 6{X,x) is continuous and strictly increasing (resp. strictly decreasing). If 
(A'a,A > 0) is SCM, then so is {9{X,X\),X > 0). 

3) Let AT be a real valued r.v. and let a : R+ —>■ Rlj_ be non-decreasing. Then (a(A) AT, A > 0) 
is SCM. 

4) Note that (SCM) implies the conditional monotonicity as defined in [18] and [5]. Indeed, 
if we take /fe = 1 for every fc € N*, we recover the conditional monotonicity hypothesis. 
The converse is not true. For example, the Gamma subordinator is conditionally monotone 
(see [5], Section 2 or [18], Section 1.4), but not strongly conditionally monotone. 

We exhibit an important class of SCM random vectors. Indeed, we prove that MTP 2 
random vectors are SCM. This result extends Theorem 3.43 of [4] to the non-Markovian 
case. 


Theorem 3.4. Every MTP 2 random vector is SCM. 

Proof. Let (Xi,--- ,X„) (n > 2) be a MTP 2 random vector, and let p : R" — >■ R+ 
denote its density. By regularisation, we may assume (without loss of generality) that p is 
continuous and strictly positive. Since MTP 2 property is invariant under permutations, it 
suffices to prove (SCM) for i = n. 

Let /fe : R —>■ Rl^, k = I,-- - ,n be continuous and strictly positive functions satisfying 
(3.1). By truncature, we may suppose (without loss of generality) that all fk are bounded. 
We shall prove by induction that, for every I G {2, • • • , n}, and every (f G Ii, 


E 

n fk{Xk) 

k^l 

Xl — Zl^ • • ‘ , Xrt — 

j ' ' 

E 

1_1 

1 1 

II 

II 



(EO 


is componentwise non-decreasing. 

Observe that if / = n, we recover (SCM) with i = n. 
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• Case 1 = 2. For every (j> G I 2 and every z = {z 2 , • • • , Zn) G R” we define: 

E [<f,(Xi,X2)fl(Xi)f2(X2)IX2 =Z2,--- ,Xn=Zn] 


/C(2,z) = 


Elfi(Xi)f2(X2)IX2=Z2, 

(f>(x,Z 2 )f(x)p(x,z)dx 

00 • 

f{x)p{x,z)dx 


1 


(3.2) 


We set 


f{y)p{y,z)dy 


F^{x) := 


/ f{y)p{y,'^)dy 

J — 00 

so that (3.2) may be written: 

/ OO 

(j){x,Z2)dF^{x) = / 4>{F~'^{u),Z 2) du, 

-OO J 0 


since F^ is continuous and strictly increasing. 

Therefore, for every u G [0,1], it suffices to show that 

z I —> F~^(u) is componentwise non-decreasing. (3.3) 

But (3.3) holds as soon as, for every a; G R., z 1 —Fz(x) is non-increasing with respect 
to each argument. Indeed, for every z = (z 2 , • • • , 2 „), z' = (^ 2 , • • • , z(j) in R”“^ such that 
z < z' (i.e. Zi < z\ for every i = 2, • • • , n), and for every u G [0,1], 

Fz' {F-\u)) < F, (F-i (u)) (since z 1 —>■ F^ix) is componentwise non-increasing) 

= u = F,, {F-\u)) . 

Since p is MTP 2 , then, for every z < z', 


p{y, z')p{x, z) > p{y, z)p{x, z') iiy>x, 


and 

p{y, z)p{x, z') > p{y, z')p{x, z) if y < x. 
Then, for every z < z' in R"“^, 


1 

F^ix) 


1-h 


< 1 + 


f{y)p{y,z)dy 

f{y)piy,z)dy 


= 1-h 


' —OO 

/*oo 


f{y)p{y,'^')p{x,z)dy 


f{y)p{y,z)p{x,z')dy 

fiy)piy,z)p{x,z')dy 


f{y)p{y,z')dy 


-= 1 - 1 - - 

fiy)piy,Z'')pix,z)dy / f{y)p{y,z')dy 

J —OO 


1 

F’z' {x) ’ 


which proves that, for every a: G R, z 1 —> F^ix) is componentwise non-increasing, and 
then (3.3) holds. 
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• Case 1 > 2. Suppose that (E;_i) holds for every function in 2/_i. Let us prove (E/) for 
every fixed (j) G Ii- 

For every z = {zi, ■ ■ ■ , Zn) G we define: 


E 


/C(l,z) = 


,xi) n fk{Xk) 

k=l 


Xl — Zlj‘ • • , X^ — Zn 


E 


n fk{Xk) 

Lfe=l 


Xl — Zl^‘ ‘ ‘ , Xn, — Zn 


/ / (j){x,xi-i,zi)W fk{xk)p{x,xi-i, 2 .)dx\ dxi-i 

J-^ \Jw-^ _^_ J 

/ / Y\ fkixk)p(x,xi-i,z)dx] dxi-i 


where x = (a;i, • • • , xi- 2 ) and dx = dxi ■ ■ ■ dxi- 2 - 

Now, set X = {Xl, • • • , Xn), and consider the functions (p : —>• K. and p : —>• 

K.^ defined respectively by: 


1-2 


4>{x, u, Zi) fk{xk)p{x, u, z)dx 


p{i 


k=l 


1-2 


fk{xk)p{x,u,z)dx 


k=l 


E 


1-1 


PiX,,--- ,Xi_,,zi) n fk{Xk) 

k=l 


Xi^i = u,X = z 


E 


n fk{Xk) 

k=l 


X /_1 = u,X = z 


and 


. 1-2 

p{u,z)= / T\ fk{xk)p{x,u,z)dx. 

fc=i 


By the induction hypothesis (E/_i), p belongs to In- 1 + 2 , and, using Point 2) of Theorem 
2.24, p satisfies (MTP 2 ). Moreover, 


/C(/,z) = 


1 , z)//_ 1 (a;/_ 1 )p(a;/_ 1 , z )1 


/ OO 

fi-i{xi-i)p{xi-i,z)dxi-i 

-OO 


Since p is MTP 2 , 
for every y GR, 


then, using the same computations as in the Case I = 2, we show that, 

pOO 

/ //_i(a;i_i)p(a;/_i,z)(ix/_i 

z I_. h. _ 

/ fi-i{xi-i)p{xi-i,z)dxi-i 


J —OO 

is non-decreasing with respect to each argument; which yields that z 1 —>■ X{l,z) is also 
non-decreasing with respect to each argument. □ 


Remark 3.5. We failed to find a SCM process which is not MTP 2 due to that SCM 
and MTP 2 properties coincide for several processes. For example, let (Xi,X 2 ) be a SCM 
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random vector which has the law P(Xi = 1 ^X 2 = j) =p{i,j) (*,J G Z). Then, (Xi,^ 2 ) is 
MTP 2 . To prove this, we may assume without loss of generality that p is strictly positive. 
Since (Xi,^ 2 ) is SCM, then, for every bounded and strictly positive function /, and for 
every a G Z, 

+00 

^/(fc)p(fc,n) 

n I —> - is non-decreasing, 

X! f{k)p{k,n) 

k— — oo 

i.e. for every n < n', every a G Z and every bounded and strictly positive function 
/ : Z —>■ R+, 


+ 00 


k'—a 


f{k) [p{k, n)p{k', n ) - p{k, n )p{k', n)] 


\k——co 


> 0 


which is equivalent to 

+00 / a— 1 \ 

E h{k) [p{k,n)p{k',n ) — p{k,n)p{k',n)]\ >0 (3-4) 

k'—a \fc= —00 / 

for every bounded and strictly positive functions g : |a, -|-oo|—>■ 1R+ and — c», a — 1] —>■ 
]R_i_. But (3.4) implies 


V fc' G |a, -l-oo|, V — 00 , a — I] —>■ M+, 

a—1 

h{k) [p{k, n)p{k', n') — p{k, n')p{k', n)] > 0 


which in turn implies 

V/c' G |a, -|-oo|, V/c g] — 00 , a — I], p{k, n)p(k', n') — p{k, n')p(k', n) > 0. 

Since a, n and n' are arbitrary integers, we deduce that p is TP 2 . As a consequence, every 
SCM process taking values in a discrete subset of R has TP 2 bidimensional marginals. In 
particular, a Markov process with discrete state space is SCM if and only if it is MTP 2 . 

Here is a direct consequence of Theorem 3.4. 

Corollary 3.6. Every stochastic process with MTP 2 finite-dimensional marginals is SCM. 
In particular, if X := {Xt,t > f) is a real valued Markov process such that X has an 
absolutely continuous and TP 2 transition kernel, then X is SCM. 

Remark 3.7. By Corollary 3.6, the processes below are SCM: 

i) processes with independent and log-concave increments, 

ii) absolute values of processes with independent, symmetric and PFqo increments, 

iii) one-dimensional diffusions with absolutely continuous transition kernel. 

iv) Gaussian random vectors with an invertible covariance matrix such that the inverse 
matrix has negative off-diagonal elements. 
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An application of multivariate total positivity to peacocks 


4 Applications of Strong conditional monotonicity to 
peacocks 

We use strong conditional monotonicity results to study some generalisations of the Carr- 
Ewald-Xiao theorem (see (CEX08)). 


4.1 Peacocks obtained by integrating with respect to a finite pos¬ 
itive measnre 

The following result was proved in [5]. 

Theorem 4.1. Let (X\,X > 0) be a real valued right-continuous process which is con¬ 
ditionally monotone in the sense that, for every n € N*, every i € {I,-"’ jn}, every 
0 < Ai < • • • < Xn, and every </> € 

z I —> E[(/)(Xai , • ■ ■ : non-decreasing. (CM) 

Suppose that, for every t > 0 and every compact K C M+, 


E 


exp 


(t sup Xa 
\ agk 


< oo and inf E[exp(tXA)] > 0. 


(4.1) 


Then, for every finite positive measure p on 1R+, 

V := “ P(^o,cock. 

This result is a generalisation of Carr-Ewald-Xiao theorem. Indeed, by making the 
change of variable s = tX in (CEX08), the Brownian scaling property yields: 

Vt > 0, N* = /' ^ (4.2) 

Jo Jo 

Observe that {B\,X > 0) is a right-continuous conditionally monotone process since it is 
a Li'^Avy process with log-concave increments. Then, by Theorem 4.1, 





>0 E 




dX, t > 0 


is a peacock. 


and we recover (CEX08) thanks to (4.2). 

In [5], further examples of conditionally monotone processes are presented. Eor example, 
the Gamma subordinator and “well-reversible" diffusions at a fixed time are conditionally 
monotone. We refer to [5] for the definition and some properties of “well-reversible" diffu¬ 
sions. Moreover, Theorem 4.1 applies to stochastic processes with MTP 2 finite-dimensional 
marginals (such as one-dimensional diffusions) since they satisfy (SCM) (which implies 
(CM)). Note that one-dimensional diffusions are not necessarily “well-reversible" at a 
fixed time. Indeed, “well-reversible" diffusions at fixed time are unique strong solutions of 
stochastic differential equations. 


4.2 Peacocks obtained by normalisation. 

Let {Vt,t > 0) be an integrable real valued process with a strictly positive mean, i.e. 
E[|Vt|] < 00 and E[Vt] > 0. Consider the process 



Vt 

E[y] 


,t > 0 
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Observe that E[A^j] = I for every t > 0. Since E[A^t] does not depend on t (which is a 
necessary condition to be a peacock), it is natural to look for processes {Vt,t > 0) for which 
{Nt,t > 0) is a peacock. Many examples of such processes are presented in [6]. Note that 

if Vt := / then (CEX08) is equivalent to: 

Jo 



Vt 

E[Vi] 


t > 0 


is a peacock. 


One may also investigate on processes (kt,t > 0) such that the centered process (Ct '■= 
Vt — E[Vt], t > 0) is a peacock. We do not treat this case here and refer to [6] and [18] for 
main results. 

We deal with processes of the forms 


Vt = exp 


q{X,Xx)^iidX) 


(Fi) 


and 


Vt^ = exp 


q{X,tXx)fi{dX) 


(F 2 ) 


Vo / 

where /i is a positive Radon measure, q : M+ x M —>■ M a continuous function such that, 
for every s > 0, : a; 1 —> q{s,x) is non-decreasing (resp. non-increasing), and where 

{X\, A > 0) is a real valued process. Our purpose is to answer the following question: 


Under which conditions is 



E[U/] 


u>0 ( 


1 ,2) a peacock? 


(Q) 


Remark that, in (Fi), the parameter t is a time parameter, while in (F 2 ), it is a dilatation 
parameter. For this reason, we call peacocks with respect to maturity peacocks of type 
and peacocks with respect to volatility peacocks of type N'^. 

For processes of type N^, a partial answer to (Q) is given in [6] when X has independent 
and log-concave increments. We extend this result to real valued processes which satisfy 
(SCM). In particular, we prove under some integrability hypotheses that processes with 
MTP 2 finite-dimensional marginals solve (Q). 


4.2.1 Peacocks with respect to maturity 

Theorem 4.2. Let {X\,X >0) be a right-continuous process which satisfies (SCM). Let 
p, be a positive Radon measure and let q : K.+ x K. —>■ K. 6e o continuous function such that, 
for every t>0: 

i) yi—> q{t,y) is non-decreasing (resp. non-inereasing), 
ii) the following integrability properties hold: 


0t := exp (/r([0, t]) sup q{s,Xs)] is integrable 

0 <s<t 


and 


At :=E 


exp /r([0,t]) inf q{s,Xs) 

0 <s<t 


Then, 


Nt := 


exp q{s,Xs) fiids) 




E 


exp^y q{s,Xs) fi{ds) 


T.^ >0 


> 0 . 


is a peacock. 


(4.3) 

(4.4) 

(4.5) 


J 
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Proof. We only consider the case where y i —> q{X, y) is non-decreasing. 

Let T > 0 be fixed. 

1) We first suppose that p. has the form: 

r 

l[o.T]C?M = ajSxi, 

i=l 

where r G N, r > 2, oi > 0 ,02 > 0,..., > 0, = /^([Oi T']), 0 < Ai < A 2 < 

Xr < T, and where 6\^ is the Dirac measure at point A^. 

We show that, 

(Nn := exp (^2aiq{Xi,X\.) — h(n)^ , n G {1, 2, • • • , r}^ is a peacock, 


where 

Note that 

with 


h{n) := logE 


exp aiq{Xi,Xx 




E[Af„ - Af„_i] = 0, for every n G {1, 2, • • • , r} 

Nn - Nn-1 = Nn-1 (e^r^q(^n,Xxn)-h{n)+h(n-l) - = Nn-1 


and 


qn{y) = anq{Xn,y) - h(n) + h(n - 1). 
Then, for every convex function ^ G C, 


- 1 


E[V;(iV„)] - E[i;{Nn-i)] > E [ij'{Nn-i)Nn-i 

K{n,XxJE[Nn-i\XxJ 


= E 


- 1 


where 


K{n, z) = 


E[lf'iNn-l)Nn-l\Xx„=z] 


E[Nn-l\Xx^ = z] 
Observe that the function (p : —>• R._|_ given by: 

/n —1 


■ ■,Xn-l) = -Ip' 


exp aiq{Xi,Xi) - h{n - 1) 




belongs to Xn-i- If, for every z € M*, we define: 

U{x) = for every x G 

then, for every n G {2, • • • , r}, 

n —1 


Nn-i=e-'^^^-^^\{h{Xx, 




(4.6) 

• < 
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n —1 

E cl)iXx„...,X^^_,)l[MXx,) X^„=z 




K{n, z) = 


Note that K(n, z) is well defined since, for every n G {1, 2, • • • , r}, 


E ]^/fc(-’fAj = E exp I y^^aiq{Xi,Xx. 


< E exp ( sup g(A, X\) ( 


< E exp j sup q{\, X\) j V I 


= E exp I Q;(r) sup q{X,X\) ) V 1 

L \ 0<A<T / 

= E[0t V 1] < oo. 


By (SCM), K{n,z) is non-decreasing with respect to z. 

Now, for every n G N*, we denote by {qn)~^ the right-continuous inverse of qn and we set: 
V{n,XxJ := K{n,XxJE[Nr,-i\Xx„] - l) . 


i) if A:a„ < (9 „)-i(0), then - I < 0 and 

V{n,XxJ >K{n, (g„)-i(0)) E[fV„_i|X aJ - l) 

ii) if A:a„ > (g„)-i(0), then - I > 0 and 


V{n,XxJ >K{n, (g„)-i(0)) E[fV„_i|X aJ 


-1 . 


As a consequence. 


E[V^(fV„)]-E,[V>(iV„-i)] 

>E[y(n,AA„)] >if (n,(g„)-i(0))E[E[N„_i|AAj - l) 

= K (n, (g„)-'(0)) E - l) 

= (n, (g„)-i(0)) E [iV„ - iV„_i] = 0; 


which shows that, for every integer r > 2, 


Nn := exp E atq{Xt, Xxi) - h{n) j ,n G {1,2, • • • ,r > is a peacock. 
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2) We consider v = l[o^T]d/r, and, for every 0 < t < T, we set: 

exp q{u,Xu)v{du) 


iv(-) = 


E 


exp^y q{u,Xu)v{du) 


Since the function A G [0,T] i —> q{X,X\) is right-continuous and bounded from above by 
sup |g(A, Xa)| which is finite a.s., there exists a sequence {vmn G N) of measures of the 

0<A<T 

form (4.6), such that, for every n G N, suppr-n C [0,T], / Un{du) = f ^{du) and, for every 

0 < t < T, 


lim exp ( / q{u, Xu)i'n{du) ] = exp ( / q{u, Xu)i'{du) ] a.s. (4-7) 

\Jo ) \Jo ) 


Moreover, for every 0 < t <T and every n G N, 

rt 


exp ^y q{u, Xu)i^n{du) 

< exp ( sup q{X,Xx) i^n{du)] 

\0<X<T Jo J 

< exp sup q{X,Xx) / Vnidu) V 1 

\ 0 <\<T Jo J 

= exp [ sup q{X,Xx) [ ) V 1 = 0t V 1 

\ 0 <\<T Jo j 


which is integrable from (4.3). By the dominated convergence theorem. 


lim E 

n—>-oo 


exp ^y q{u,Xu)vn{du^ =E exp ^y q{u,Xu)v{du) 


Using (4.7) and (4.8), we obtain: 


lim Nf"’ = Nf’ a.s., for every 0 < t < T. 

n—¥oo 


But, we proved in 1) that: 

0 < t < is a peacock for every n G N, 
i.e., for every 0 < s < t < T and every i/j € C: 


E 




< Er 


Besides, 


sup sup 

0<t<T n>0 


N, 




< 


(Brp V 1 

A 1 


(4.8) 


(4.9) 


(4.10) 


(4.11) 


(4.12) 


which is integrable from (4.3) and (4.4). Using the dominated convergence theorem, we 
pass to the limit in (4.11) as n —>■ oo and deduce that < t < T) is a peacock for 

every T > 0. □ 
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Now, we prove a version of Theorem 4.2 for a squared Bessel of dimension 0 (denoted 
BESQ°). Note that the transition function of a BESQ° is not absolutely continuous with 
respect to Lebesgue measure. Then Theorem 3.4 does not apply. In particular, the finite¬ 
dimensional marginals of a BESQ° do not satisfy (SCM). Nevertheless, a limit theorem 
due to Teller [14] for critical Galton-Watson branching processes allows to exhibit peacocks 
of type (4.5). 

Example 4.3. (A version of Theorem 4.2 for a BESQ°). 

For every k G W, let Z'^ := (Z^ , n G N) denote a Galton-Watson branching process 
starting with k individuals, and which has a geometric reproduction law v of parameter 
1 . 

2 ’ 

v{i) = 2“*“^, for every i G N. 

For every fc G N*, is an homogeneous Markov chain with values in N, and its transition 
probability matrix Q is given by: 

Vi£N, 4.(0,i) = {j 

and 

V(z,j) G N* X N, = (^ ^ ^ 2-(*+^). 

We consider the family G N) of transition functions defined on N x N by: 




1 if I = j 
0 otherwise 


and 

Vn > 1, j) = X! 

nSN 

Since the function (f, j) '—^ j ^ TP 2 , we deduce from (2.3) that is TP 2 
for every n G N. 

For every A > 0 and every fc G N*, we set: 



[feA]> 


where [•] denotes the floor function. Then, (lAf, A > O) is a Markov process with values in 
—N, and its transition function {PQ^rj,0 < C < ''l) is given by: 

fx 

'ix,yG PQ^n{x,y) = 


Observe that is TP 2 for every C, and ij. By Theorem 3.4, for every k G N*, (lAf, A > O) 
is SGM. 

Let q : R+ x R —>■ R be bounded, continuous and such that, for every A > 0, a; 1 —>■ g(A, x) is 
non-decreasing (resp. non-increasing). It follows from Theorem 4.2 that, for every > 0, 
i G N*, and for every strictly increasing sequence {Xi,i> 1) in R^, 


exp E aiq{\i,Yy 


N: := 




E 


exp i'^aiq{Xi,Yj^J 


■,n G N 


is a peacock. 


(4.13) 
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A result due to Feller [14] states that, as k tends to oo, {Y^, A > O) converges in distribution 
to (Y^, A > 0), which is the unique strong solution of: 


dZ\ = ^2Z\dB\, Zo = 1, 

where {B\, A > 0) is a standard Brownian motion. In particular, (lAf, A > O) converges in 
sense of hnite distributions to A > 0). Then, (4.13) yields: 


exp ( ^a*( 7 (Ai,r;(^) 


:= 




E 


exp 


'^a,q{Xi,Y) 


OO \ 
Ai ) 




■,n G N 


is a peacock. 


(4.14) 


As a consequense, we obtain the following result: 

Corollary 4.4. Let (Yt,t > 0) be a BESCf issued from 1, and let q : K_|_ x R+ — >■ K &e a 
continuous and bounded function such that, for every A > 0, y i —> q(\, y) is non-decreasing 
(resp. non-increasing). Then, for every positive Radon measure pL on K.+ , 


/ 


Nt 


V 


exp^y q{s,Ys)p.{ds) 


E 


exp (y q{s,Ys)fi{ds) 


0 


is a peacock. 


4.2.2 Peacocks with respect to volatility 

Theorem 4.5. Let {X\, A > 0) a right-continuous process which is SCM. Let q : IR+ x K. — >■ 
R 6e a continuous function such that, for every A > 0, a; i —> q{X,x) is of class, and let 
pL denote a positive Radon measure on M_|_ satisfying: 


Vt > 0, E 


exp 



q{X,tXx)piidX) 


< oo. 


We suppose that: 


i) 


for every A > 0, the functions x i —> q{X, x) and x 
(resp. non-increasing). 


a;|^(A,a;) 


are non-decreasing 


ii) for every t,X> 0, there exists a = aft, A) > 1 such that: 


E 


{x,tx,) 


< oo. 


Hi) for every t,j5 > 0, and for every compact K C M+, 

:= exp ( P sup q{X,tX\) ) is integrable. 


and 




exp ( /3 m^q{X,tXx) 


> 0 . 


(4.15) 


(4.16) 


(4.17) 
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Then, 


■= 


exp q{X,tXx)^i{dX) 




V 


E 


exp^y q{X,tX\)fj,{dX) 


,t>0 


is a peacock. 


(4.18) 


) 


Proof. We shall suppose without loss of generality that, for every A > 0, the functions 
X I —> q{X,x) and x i —> x-^{X,x) are non-decreasing. 

1) We first treat the case where /r is of the form 

n 

i=l 


where m G N*, Oi > 0, • • • , a„ > 0, 0 < Ai < • • • < A^, and where denote the Dirac 
measure at point A. Precisely, we show that 


Nt := exp E aiq{Xi,tX\f) — h{t) j , t > 0 J is a peacock. 




with 


h{t) = logE 


exp E aiqifXi, tX \^) 




dq 


We set /i := X) Since the functions x i —> q{X,x) and x i —> x—{X,x) are non- 

i=l 

decreasing, then, for every 0 < b < c, and every t € [b, c ], 




and for every i € {1, • • • , m}, 


\XxA^{K,tXx,) < ^\XxA^iXi,cX^,). 


We deduce from (4.15), (4.16), (4.19) and (4.20) that, for every 0 < & < c. 


E 


I ^ dq ( 

sup < y^ai\X),^\-^{Xi,tXxi)exp I y^akq{Xk,tXx,^) 

_te[b,c] X \k=i / 


< oo. 


exp aiq{Xi,tXx.) < exp p, sup g(Ai,0) -f exp p, sup q{Xi,cXx.) , 
\~! / \ J y iG{l,--,m} J 

(4.19) 

(4.20) 


(4.21) 


Consequently, h is continuous on [0, -l-oo[, differentiable on ]0, -l-oo[, and for every t > 0, 


h'(t)e'‘W =^a,E 




Xx^^{Xi,tXx.)exp I y^afcg(Afc,tXAj 


‘ dx 


Kk=l 


i.e. 


Now, define 


h'{t) = ai E 

i=l 

hxi{t) = E 


dq 

NtXx,^{Xi,tXx,) 


NtXxAx,,tXx,) 

ox 


(4.22) 

(4.23) 
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so that 


h'{t) = y^^Qjhx^ (t). 


i=l 


Since E[Nt] = 1, then, for every t > 0 and i G {1, • • • , n}, (4.23) yields: 


E 


do — 

Nt[Xx,-^{X,,tXx,)-hx,{t) 


= 0 . 


(4.24) 


(4.25) 


On the other hand, if is a convex function in C, then (4.16), (4.22) and (4.24) imply 


d 


-E[i;{Nt)] = J2a.E 


{ Xx,^{\,tXx,) - hxM 


Thus, it remains to prove that, for every i G {1, • • • , m}. 


dq 


( Xx,-^{Xi,tXx,) - hxM 


> 0 . 


(4.26) 


A, :=E 
Observe that the function 

(j) : {xi ) I—■*/'' ^exp akq{\k, txk) - h{t)^ ^ 

belongs to Moreover, if, for every k G {I,-"’ we set fk{x) = ex.p{akq{Xk,tx)), 

then 


iVt=n 




Therefore, by setting 

Ki(m,z) := 

for every 2 ; G K., and every t G {1, • • • , m}, we obtain 


E 

< 

1 _ 1 

m 

■■ ,xxj n fk{xx,) 

k=l 

1 -1 

II 

E 

m 

n fkixx,) 

.k^l 

Xx,=z 



A, =E 


K,im,XxMNt\Xx,] ( Xx,^iXi,tXx,) - hxM 


dq 


By hypothesis i), qr. : x 1 —> x—(Q, tx) — hr. (t) is continuous and non-decreasing; let m 

ox 

denote its right-continuous inverse. Since (Xa, A > 0) is SCM, the function z 1 — > Ki{'m, z) 
is non-decreasing, and we deduce from (4.25) that: 


A, >iG,(m,5);-/(0))E 
= AT* ( 771 , 57 /( 0 )) E 


E[Aft|XAj (^Xx,^{K,tXx,) - hxM 

NtUx,^{Xi,tXx,)-hx,{t)\ \ =0. 


Thus, (A"t,t > 0) is a peacock. 

2) If fx has a compact support contained in a compact interval of K+, then, following the 
same lines as Point 2) in the proof of Theorem 4.2, we prove that ^0^ is a 

peacock. 
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3) In the general case, we consider the sequence (/r„((iA) := l[o,„]/r(dA), n G N). Let tp G C. 
By Point 2) above, ^ peacock for every n. Then 


VO < s < t, E 




(4.27) 


Moreover, it follows from Theorem 4.2 that, for every t > 0, 


( 




exp(y q{C,tX^)iJ,{dC) 


\ 


V 


E 


exp(y q{C,tXi^)fi{dC) 


,n > 0 


is a peacock. 


J 


in other terms, the sequence ^E ,n > 0^ is non-decreasing and bounded from 

above by E . Therefore, letting n tends to oo in (4.27), we obtain: 

VO < s < t, E < E hiM'"'’)] 


which proves that ,t > 0^ is a peacock. □ 

We end with some examples. 

Example 4.6. Let /i be a positive Radon measure. Let (AT^, A > 0) be a right-continuous 
process having MTP 2 finite-dimensional marginals, and such that, for every /3 > 0, 


and 


E 


E 


exp j3 sup X\ 

\ o<c<i , 


exp ( (3 inf Xr 

' 0<A<1 


< 00 


> 0 . 


We consider the function q : [0,1] x R —>■ M defined by: 


V (A, x) G [0,1] X R, q{\, x) = 2x + \/l + X + x"^. 


The following inequalities are immediate. 

V (A, x) G [0,1] X R, . 

Then, using (4.28) and (4.30), we have: 


(4.28) 


(4.29) 


(4.30) 


Vt > 0, 


E 


exp 


m([ 0,1]) sup q{X,tXx) 

0<A<1 


< OO. 


Moreover, (4.28), (4.29) and 
fulhlled. 

On the other hand, x 1 —> 
Therefore, by Theorem 4.5, 


(4.30) ensure that conditions (4.15)-(4.17) of Theorem 4.2 are 


q(\, x) and x 




are non-decreasing functions. 


( 

Nt 

V 


exp ^y q{X,tXx)fi{dX)j 


E 


exp ^y q{X,tXx)fiidX) 


\ 

■T,t> 0 


is a peacock. 
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Example 4.7. Let X := (Xu,u > 0) denote a right-continuous process which has MTP 2 
finite-dimensional marginals. We assume that X enjoys the scale property of order 7 > 0, 
i.e. 

Vt > 0 , {Xtu,u > 0) {t'^Xu,u > 0). 

Let g : R —>■ R be a non-decreasing C^-function such that: 

i) the function x 1 —> xq'{x) is non-decreasing. 

ii) For every /3,t > 0 and a; G M: 


E exp P sup q{f^Xu) ) < 00 et E exp (/3 inf q{t'^Xu) ) > 0. 


iii) For every t,u > 0, there exists a = a{t, u) > 1 such that 

E[\xx{qTit^x^)]<^- 

Then, after the change of variable s = tu, we deduce from Theorem 4.5 that 


Nt := 


qi^s)ds^ 
E exp Q J q{Xs)ds 


t > 0 is a peacock. 


The purpose of the next example is to exihibit peacocks of type using processes 
which are not MTP 2 . 

Example 4.8. Let {Xi,X 2 ) be a random vector with values in {1,2,3}^ and which has 
the law P(A'i = i, X 2 = j) = Pij, where 

1/3 3 1\ 

P= — \ 3 2 2 . 

\ 1 2 3 / 

Observe that P is a symmetric matrix which is not TP 2 since 


But, (Xi,X 2 ) satisfy a SCM type condition (weaker than SCM). Precisely, for every 
componentwise non-decreasing function p : {1,2,3}^ —>■ R and every non-decreasing and 
strictly positive functions fi, f 2 ■ {1,2,3} —>■ R+, the maps Ki, K 2 : (1, 2,3} —>■ R defined 
by 

r. ■ . r. o o. , E[0 (Xi,X2)/i(Xi)/2(X2)]|Xi = i] 

K...e{l,2,3)^ -E|;,,X,)/,(x,)|X,=i]- 


K2:i€ {1,2,3} 


E[(j){XiX2)h{Xi)h{X2)]\X2 = i] 


' ‘ E[h{X^)f2{X2)\X2=i] 

are non-decreasing. 

Since P is symmetric, it suffices to prove that ATi is non-decreasing. Observe that, for 
every i = 1,2, 3, 

I] (l>{i,k)f2{k)Pik 

KX) = -■ 

E h{k)P,u 

fc=i 


25 



ANTOINE MARIE BOGSO 


Moreover, to show that Ki is non-decreasing, we may restrict ourselves to the functions 
(j)(i, •) = I|o,-i-oo[ (a = 2, 3). Precisely, it is sufficient to see that P satisfies 


1 

Wi 


k=2 


and 


P^ 


13 


1 

Wl 

< 


k=2 


P 


23 


< 


^31 k=2 

(4.31) 

P 33 

(4.32) 


4 f2{k)Plk 4 f2{k)P2k 4 f2{k)P3k 

fc=l 

To obtain (4.31) and (4.32), one may remark that, for every i € {1,2,3} and a € {2,3}, 


4 f2ik)Pk 

k—a _ ^ 

E f2{k)Ptk 4 f2{k)Pik 

fe=l _|_ k=l _ 

4 f2ik)Pk 

k—a 

Now, since 

^[3/2(2) + /2(3)] < ^[/2(2) + /2(3)] < 2/2(2) + 3/2(3) (4.33) 

and 

1 ^ 2 ^ 3 

3/2(1) +3/2(2) - 3/2(1) +2/2(2) - / 2 ( 1 ) + 2/2(2)’ 
we deduce that Ki is non-decreasing. Note that the first inequality in (4.33) holds since 
/2 is non-decreasing. 

Similarly, one may prove that if (Yi, Y2) is a random vector taking values in {1,2,3}^ and 
whose law is given by P(Yi = z, Y2 = j) = Pij, with 


P* 


1 


3 2 1 \ 

2 2 3 
13 3/ 


then, for every componentwise non-decreasing function (p : {1,2,3}^ —>■ R and for every 
non-increasing and strictly positive funetions 51,52 : {!) 2,3} —>■ R, 


K* 


E[0(yi,y2)gi(Yi)g2(Y2)|Yi =z] 

E[gi(yi)g2(Y2)|Yi = i] 


and 


are non-decreasing. 


E[/.(yi,y2)gi(Yi)g2(Y2)|Y2 = i] 

E[gi(Yi)g2(Y2)|Y2 = i] 


Corollary 4.9. If qi,q 2 : R+ —>■ R are two non-decreasing class functions such that 

X I — > xq'i{x) and x 1 — > xq2{x) are also non-decreasing, then 


and 


are peacocks. 



exp{qi{tXi) -\- q2{tX2)) ^ 

E[exp(q'i(tYi)-I-92(1^2))] ’ “ 



exp(-gi(tYi) - g2(tY2)) ^ „ 

E[exp{-qi{tYi) - q 2 {tY 2 ))]' ~ 


26 



An application of multivariate total positivity to peacocks 


Acknowledgements: We are grateful to Professors Francis Hirsch, Bernard Roynette 
and Marc Yor for fruitfull suggestions and comments. We also thank the referee for his 
comments and technical remarks. 


References 

[1] M. Y. An. Log-concave probability distributions: Theory and statistical testing. SSRN, 
pages i-29, May 1997. 

[2] R. H. Berk. Some monotonicity properties of symmetric Polya densities and their 
exponential families. Z. Wahrscheinlichkeitstheorie und Verw. Gebiete 42: 303-307, 
1978. 

[3] D. Baker, C. Donati-Martin, and M. Yor. A sequence of Albin type continuous martin¬ 
gales, with Brownian marginals and scaling. In SA@minaire de ProbabilitA@s XLIII, 
44 I- 449 , Lecture Notes in Maths., 2006, Springer, Berlin, 2010. 

[4] A. M. Bogso. Etude de peacocks sous des hypothA "ses de monotonie conditionnelle et 
de positivitA@ totale ThA”se de I’UniversitA© de Lorraine 2012. 

[5] A. M. Bogso, C. Profeta and B. Roynette. Some examples of peacocks in a Markovian 
set-up. In SA@minaire de ProbabilitA@s, XLIV, Lecture notes in Math. (2046) 281- 
315, Springer, Berlin (2012). 

[6] A. M. Bogso, C. Profeta and B. Roynette. Peacocks obtained by normalisation, strong 
and very strong peacocks. In SA@minaire de ProbabilitA@s, XLIV, Lecture notes in 
Math. (2046) 28I-3I5, Springer, Berlin (2012). 

[7] D. Baker and M. Yor. A Brownian sheet martingale with the same marginals as the 
arithmetic average of geometric Brownian motion. Elect. J. Prob., 14(52)1532-1540, 
2009. 

[8] D. Baker and M. Yor. On martingales with given marginals and the scaling property. In 
SA@minaire de ProbabilitA@s XLIII, 437-439, Lecture Notes in Maths., 2006 Springer, 
Berlin, 2010. 

[9] P. Carr, C.-O. Ewald, and Y. Xiao. On the qualitative effect of volatility and duration 
on prices of Asian options. Finance Research Letters, 5(3):162-171, September 2008. 

[10] H. Daduna and R. Szekli. A queueing theoretical proof of increasing property of Polya 
frequency functions. Statist. Probab. Lett., 26(3): 233-242, 1996. 

[11] A. Edrei. On the generating function of a doubly infinite, totally positive sequence. 
Trans. Amer. Math. Soc. 74(3): 367-383, 1953. 

[12] B. Efron. Increasing properties of Polya frequency functions. Ann. Math. Statist., 
36:272-279, 1965. 

[13] M. Emery and M. Yor. A parallel between Brownian bridges and gamma bridges. 
Publ. Res. Inst. Math. 5'ci.40(3):669-688, 2004. 

[14] W. Feller. Diffusions processes in genetics. Proc. Second Berkeley Symp. Math. Statist. 
Prob., University of California Press, Berkeley, pp. 227-246, 1951. 

[15] P. J. Fitzsimmons, J. W. Pitman, and M. Yor. Markovian bridges: construction. 
Palm interpretation, and splicing. Seminar on Stochastic Processes, 1992 (Seattle, WA, 
1992), IOI-I 34 , Progr. Probab., 33, Birkhauser Boston, Boston, MA, 1993. 


27 



ANTOINE MARIE BOGSO 


[16] R. D. Gupta and D. St. P. Richards. Multivariate Liouville distributions. J. Multi¬ 
variate Anal. 23:233-256, 1987. 

[17] K. Hamza and F. C. Klebaner. A family of non-Gaussian martingales with Gaussian 
marginals. J. Appl. Stock. Anal, pages Art. ID 92723, 19, 2007. 

[18] F. Hirsch, C. Profeta, B. Roynette, and M. Yor. Peacocks and associated martingales. 
Bocconi-Springer, vol 3, 2011. 

[19] F. Hirsch and B. Roynette. A new proof of Kellerer Theorem ESAIM: PS, 16:48-60, 
2012 . 

[20] F. Hirsch, B. Roynette, and M. Yor. From an ltd type calculus for Gaussian processes 
to integrals of log-normal processes increasing in the convex order. J. Math. Soc. Japan, 
63(3):887-917, 2011. 

[21] F. Hirsch, B. Roynette, and M. Yor. Unifying constructions of martingales associated 
with processes increasing in the convex order , via Levy and Sato sheets. Expositiones 
Mathematicae, 4, 299-324, 2010. 

[22] F. Hirsch, B. Roynette, and M. Yor. Applying Ito’s motto: “look at the infinite 
dimensional picture” by constructing sheets to obtain processes increasing in the convex 
order. Periodica Mathematica Hungarica, 61(1-2):195-211, 2010. 

[23] S. Karlin. Total positivity, absorption probabilities and applications, Trans. Amer. 
Math. Soc. Vol. Ill, pp. 33-107, 1964. 

[24] S. Karlin. Total positivity. Stanford University Press, 1967. 

[25] S. Karlin and J. L. McGregor. Goincidence probabilities. Pacific J. Math. Vol. 9, pp. 
1141-1165, 1959. 

[26] S. Karlin and J. L. McGregor. Classical diffusion processes and total positivity, J. 
Math. Anal. Appl. Vol. 1, pp. 163-183, 1960. 

[27] S. Karlin and H. M. Taylor. A second course in Stochastic processes. New York: 
Academic Press, 1981. 

[28] H. G. Kellerer. Markov-Komposition und eine Anwendung auf Martingale. Math. 
Ann., 198:99-122, 1972. 

[29] S. Karlin and Y. Rinot. Classes of orderings of measures and related correlation in¬ 
equalities. I. Multivariate totally positive distributions. J. Multivariate Anal. 10(4):467- 
498, 1980. 

[30] A. Muller and M. Scarsini. Stochastic comparison of random vectors with a common 
copula. Mathematics of Operations Research, 26(4):723-740, 2001. 

[31] G. Pages. Functional co-monotony of processes with an application to peacocks and 
barrier options. In Seminaire de Probabilites XLV, Lecture Notes in Maths. Springer. 
To appear 2013. 

[32] M. Rothschild and J. E. Stiglitz. Increasing risk. I. A definition. J. Econom. Theory, 
2:225-243, 1970. 

[33] D. Revuz and M. Yor. Continuous martingales and Brownian motion, volume 293 
of Grundlehren der Mathematischen Wissenschaften [Fundamental Principles of Math¬ 
ematical Sciences]. Springer-Verlag, Berlin, third edition, 1999. 


28 



An application of multivariate total positivity to peacocks 


[34] T. K. Sarkar. Some lower bounds of reliability. Technical Report, No 124, Dept, of 
Operations Research and Statistics, Stanford University, 1969. 

[35] I.J. Schoenberg. On Polya frequency functions 1. The totally positive functions and 
their Laplace transforms. J. Analyse Math. 1, 331-374, 1951. 

[36] M. Shaked and J.G. Shanthikumar. Stochastic orders and their applications. Proba¬ 
bility and Mathematical Statistics. Academic Press, Boston, 1994. 


29 



